Electron Pair Resonance in the Coulomb Blockade 
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We study many-body corrections to the cotunneling current via a localized state with energy 
at large bias voltages V. We show that the transfer of electron pairs, enabled by the Coulomb 
repulsion in the localized level, results in ionization resonance peaks in the third derivative of the 
current with respect to V, centered at eV — ±2ed/3. Our results predict the existence of previously 
unnoticed structure within Coulomb-blockade diamonds. 

PACS numbers: 73.23.-b, 73.23.Hk, 73.63.Kv 



Introduction. — Current flow through a single localized 
state (LS) coupled to metallic leads is a paradigm of 
quantum transport through nanostructures, with appli- 
cations to many systems such as impurities embedded 
in tunnel barriers, quantum dots, single-molecule junc- 
tions, or carbon nanotubes, see, e.g., Refs. P, [3, H, 0, Q • 
Despite its simplicity, this system exhibits a wide range 
of transport behaviors, including resonant and sequential 
tunneling, cotunneling, and the Kondo effect. 

All of these regimes are captured remarkably well by a 
simple extension of the Anderson impurity model 
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FIG. 1: (Color online) Schematic stability diagram of a 
single- level quantum dot. The thick lines within the shaded 
Coulomb blockaded region are characterized by a resonance 



■ [tadl^ci^aa + ^a'^kaa'^'A ' peak in d^I/dV^ due to opening of two-electron ionization. 



which describes tunneling of amplitude ta between the 
spin-degenerate LS of energy (with creation operator 
'^1=1 I ^-iid number operator Ua = rfj-dcr) and two leads 
a = L,R (with dispersion Ck and creation operator cj^.^^). 
For large on-site Coulomb repulsion U, double occupa- 
tion of the LS is suppressed and the nature of transport 
depends on both Cd (tunable by a gate voltage V^) and the 
bias voltage V. Within the shaded areas of the stability 
diagram in Fig.[Tl the average occupation Ud = n| -|-n| of 
the LS is close to integer and current flow is suppressed 
by the Coulomb blockade. In contrast, current can flow 
by sequential tunneling processes outside the shaded ar- 
eas, where the average occupation of the dot is no longer 
integer. This picture of the Coulomb blockade has been 
confirmed in numerous experiments performed on various 
systems 

It is the main point of this paper that even the minimal 
model of Eq. ([1]) predicts additional structure within the 
Coulomb-blockaded region, emerging from two- electron 
ionization of the LS at large biases. This ionization pro- 
cess is an effect of many-body correlations, enabled by 
the on-site Coulomb repulsion, which is much more ro- 
bust than the Kondo correlations emerging in the Kondo 



For eV below the threshold voltage eVc = 2ed/3, tunneling of 
electron pairs is a precursor effect to two-electron ionization. 

valley Ud = 1 at low temperatures and small voltages. 
Indeed, the fine structure due to the two-electron ioniza- 
tion process exists in both Kondo and non-Kondo valleys, 
as illustrated in Fig. [T] 

The two-electron ionization requires biases beyond a 
threshold voltage Vc, indicated by the thick black lines 
in Fig. [T] Below the threshold voltage, correlated two- 
electron transfers between the two leads constitute a pre- 
cursor effect to two-electron ionization. While the limit 
of the Coulomb blockaded region is characterized by a 
resonance peak in the differential conductance dl/dV , 
we find that the onset of two-electron ionization at Vc is 
accompanied by a peak in d^I/dV^. Interestingly, the 
difference between both resonance phenomena emerges 
solely from familiar Fermi liquid phase space factors 
which appear in the two-electron ionization rate. One 
important implication of this analogy is that the onset 
of two-electron ionization is accompanied by anomalous 
temperature sensitivity, even when eVc 3> T, as is famil- 
iar for the boundary of the Coulomb blockaded region. 

Most of our conclusions carry over to many-level quan- 
tum dots ("metallic dots") where the stability diagram 
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FIG. 2: (Color online) Stability diagram of a metallic quan- 
tum dot with vanishing level spacing. The threshold lines 
for two-electron ionization can be determined by examining 
the electrostatic energy of the circuit. Graphically they are 
obtained {e.g., for the central diamond) by rescaling the di- 
agonal AD by a factor 3, so that AB — BC — CD. Inset: 
Equivalent electric circuit for a metallic island. 

exhibits a sequence of Coulomb diamonds, reflecting the 
step- wise population of the dot with increasing gate volt- 
age. This is depicted in Fig.[2]where we include the effects 
of asymmetric capacitances between dot and electrodes. 

In the remainder of the paper, we quantify the behavior 
of the current near the two-particle threshold. 

Two- electron ionization. — Ionization by means of sin- 
gle particle tunneling becomes energetically allowed when 
the source chemical potential eV/2 is aligned with the LS, 
i.e., at eV = i2ed. In contrast, the two-particle ioniza- 
tion process, responsible for the predicted boundaries in 
the stability diagram, is depicted in Fig. [3l^a). At finite 
bias, an electron tunneling between the leads can suffer 
an energy loss up to eV. Due to the on-site Coulomb 
repulsion, this energy loss can be transferred to a second 
electron from the source electrode, exciting it to energies 
up to 3eI^/2. Specifically, the second electron can pop- 
ulate (and thus ionize) the LS once its maximal energy 
exceeds e^, i.e., for biases exceeding the threshold voltage 



eVc = 2ed/3. 



(2) 



The predicted lines in the stability diagram originat- 
ing from the onset of two-electron ionization occur for 
V = ±Vc. Thus, they are located within the Coulomb 
blockaded region which extends up to eV = ±2ed- 

Microscopically, the two-electron ionization process 
proceeds as follows, cf. Fig. [3ja): (i) An electron with 
energy ei from the source electrode (L) enters the LS 
and (ii) tunnels into the state Ei of the drain (R). In the 
same process, (iii) a second electron with opposite spin 
and energy £2 tunnels from the source into the LS. The 
amplitudes of the steps (i) and (iii) are proportional to 
ii, while the amplitude of step (ii) is proportional to 
Thus, the resulting amplitude of two-electron ionization 
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(ed-ei)(£;i-ei)- 



(3) 



Following standard perturbation theory, the energy de- 
nominators are given by the difference between the in- 
termediate and initial energies. In Eq. ([3]), we assumed 
a large on-site Coulomb repulsion U so that there is no 
contribution from virtual states with double occupation 
of the LS. If these states were included, the correspond- 
ing terms would exactly cancel the amplitude Eq. ^ 
in the limit of vanishing U . This makes it manifest that 
two-electron ionization is enabled by the on-site Coulomb 
interaction. 

Based on Eq. ([3]) and energy conservation, the two- 
electron ionization rate per spin, at T = 0, is 
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■S{ei+e2-ed- El), (4) 



where 1. Here = 27r|ti,pi/ and F^, = 27r|t/jpi/ are 
the partial widths of the LS due to escape to source and 
drain, respectively, and v denotes the density of states in 
the leads. Performing the integration over 62, we obtain 



~ (27r)2 



v/2 

dei 



dEi- 



-eV/2 



{eV/2 + £1 - - ed) 
[ed-eiYiEi-eif 



(5) 

where 9{x) is the step function. Since —eV/2 < Ei and 
ei < eV/2, the argument of the 9{x) function is negative 
for ed > 3eV/2, i.e., for eV < eVc- In contrast, for 
< V — Vc Vc, the integration regions for ei and Ei 
are restricted to ed — eV<ei< eV/2 and —eV/2 < Ei < 
eV — ed, respectively. Since both regions are narrow, we 
find the threshold behavior 
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(6) 



of the two-electron ionization rate Fjon- 

It is crucial that energy exchange between electrons 
in the leads does not require direct interaction between 
them. Instead, this process is enabled by the finite 
Coulomb repulsion in the LS alone. In this regard, the 
underlying physics of two-particle ionization is similar 
to that of energy exchange between electrons in a bulk 
metal, facilitated by a magnetic impurity Indeed, it 
is the non-zero on-site Coulomb repulsion U that ulti- 
mately generates the magnetic impurity Curiously, 
similar many-body processes can also be enabled by the 
pairing interaction in devices consisting of two Josephson 
junctions in series, where they lead to subgap structure 
in the current [8|. 
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FIG. 3: (Color online) Schematic rendering of (a) two-electron 
ionization of the LS, in which one electron tunnels from source 
to drain while the other jumps into the LS, and (b) pair 
tunneling, in which two electrons tunnel between source and 
drain. Both processes are enabled by the on-site Coulomb 
repulsion suppressing double occupation of the LS. 

After entering the empty LS with rate Fion by two- 
electron ionization, the electron rapidly escapes into 
source or drain electrode by single-electron tunneling. 
These depopulation processes occur with rates F^ and 
Fi^, respectively. Thus, the average occupation of the LS 
is governed by the rate equation 



2Fio„(l - - nx) = (F^ + F^Jn^. 



(7) 



Here, the factor 2 accounts for spin. Eq. ([7]) yields 
rid — 2Fion/ (r^ + Lij). Since the net charge transfer is 
2e (e) when the electron tunnels out to the drain (source) 
electrode, the "two-electron ionization" current 2{V) be- 
tween the leads becomes 



X(y) = e(2Ffl + Fi)rirf«2eFi, 



2F, 



Tr 



(8) 



Due to Fion, the ionization current I{V) also exhibits the 
threshold behavior I{V) x {V - Vcf9{V - K)- 

Clearly, the ionization current, Eqs. ([8]) and ([6]), con- 
stitutes but a small fraction ^ Tl{V — Vc)^ /V^ of the 
cotunneling current ^ '"^T" V. Thus, it is an impor- 
tant question how the threshold anomaly Eq. ([S]) can be 
distinguished from the background cotunneling current. 
Eq. ([8]) predicts that two-electron ionization induces a 
jump in (PX/dV^ located a.t V = Vc- We now turn to a 
more careful analysis of this jump, focusing first on the 
two-electron current below threshold, before deriving a 
general interpolation formula. 

Two-electron current below threshold. — For voltages 
below the threshold, V < Vc, ionization of the LS is 
no longer possible by two-electron processes. But two- 
electron processes can still excite electrons in the leads 
to just below the energy of the LS. We will now show 
that this constitutes a precursor effect to two-electron 
ionization which contributes a logarithmically singular 
threshold dependence to the differential conductance. 

For large on-site Coulomb repulsion U, the two- 
electron process below threshold proceeds microscopi- 
cally as follows, cf. Fig.[3l[b): (i) A spin- up electron from 
lead ai with energy ei enters the LS; (ii) the electron tun- 
nels out to state Ei in lead a'l ; (iii) a spin-down electron 



from lead a2 with energy £2 enters the LS and (iv) leaves 
into state E2 in lead The corresponding amplitude is 



^ei-»-El 

£2 — 
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[ed - ei){Ei - ei){ed - E2) 



+ (1-2), (9) 



where the second term accounts for the four-step pro- 
cess described above with the interchanged order (iii) 1— > 
(iv)H^ (i) ^ (ii). This results in a scattering rate 

C^ah27r;.yn[de,dii;,/(e, - - f{E, - /i^'^)}] 

(10) 



1=1 

X \Al-Z%\?5{e,+e2-E,-E2) 



Here /(e) = [e'/^ + l]"i and ^^/^ = ±eV/2. The result- 
ing two-electron tunneling current contains two contri- 



2eFi: 



corre- 



butions, I = /(i'^) + where /(^'^^ 

sponds to two-electron transfer between the leads, while 
j(ie) ^ e^„(F^i:^ -I- F^^'^) accounts for one-particle 
transfer between the leads, accompanied by the creation 
of a particle-hole excitation in one lead. 

The crucial observation is that X is singular as V ap- 
proaches Vc from below. The singularity arises from the 
domain Ei ~ ~eV /2, E2 ~ ieV /2, ei ~ €2 ~ eV /2. To 
see this, we first note that in this domain, the amplitude 

(IvlyfeA^ ■ Using the 
Golden Rule Eq. (fTO]) . and performing the integrals over 
ei, e2, and Ei, we obtain for eVc S> T 



A%Zt simplifies, A^Zt 



2eTlTR{TR + \TL) 
h {2TiY{eVcY 



dE: 



l-f{E2 + eV/2) 
{ed - E2Y 



xf{E2 - 3ey/2j y{TiTf + {E2 - 3eV/2f^ . (11) 

Since both E2 and Cd in the denominator of Eq. (fTT|) 
are close to 3ey/2, the remaining integration yields the 
singular contribution 



dl 2e2 3F2 F;^(F^ + iF^) 



dV 



(27r)2(eK)^ 



In 



eVc 



max{eVc - eV, T} 



(12) 



to the differential conductance. The logarithmic singular- 
ity in the two-electron tunneling current at Vc signals the 
opening of the two-particle ionization channel in Eq. ([5]) 
which is lower order in the tunneling amplitudes and in- 
volves real occupation of the LS. 

The appearance of T-dependence in Eq. at eVc ^ 
T resembles the behavior of the conductance near the 
onset of sequential tunneling at eV = ±2ed, cf. Fig. 1. 
In fact, we find that the analogy between the onset of 
sequential tunneling at eV — ±2ed and the onset of two- 
electron ionization at V = ±Vc goes much further. The 
lines eV = i:2ed in the stability diagram separate trans- 
port regimes with real occupation (sequential tunneling) 
and virtual occupation (cotunneling) of the LS. Similarly, 
the lines V = ±Vc separate regimes with real occupation 
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(two-electron ionization) and virtual occupation (pair- 
tunneling) of the LS. We now explore this analogy on a 
quantitative level. 

Correspondence of one-electron and two- electron ioni- 
zation. — We start by noting that Eqs. ([8|) and yield 
d^I/dV^ cx 0{V - Vc) and d^I/dV^ oc 1/(14 - F), which 
are the familiar voltage dependencies of the sequential- 
tunneling and cotunneling currents, respectively, pro- 
vided we make the replacement eVc 2ed- This suggests 
that the currents near the onsets of sequential tunneling 
and two-electron ionization are related to one another 
more generally by two voltage derivatives. To establish 
this relation, although approximately, we incorporate the 
lifetime broadening T = Tl+Tr of the LS into Eq. pT|) . 
and cast it into the form 



2erlrR{TR 



^Fl) 



de{f(e-eV)-f{e + eV)} 



[e-{ea-eV/2)r + {r/2y 



(27r)2(eK)4 
(nTf + {e- eVf 



(13) 



where e = E2 — eV/2. The finite lifetime provides a 
physical cutoff of the singularity in Eq. (H^). Most im- 
portantly, Eq. p3p captures processes involving both vir- 
tual and real occupations of the LS, i.e., it describes the 
two-electron resonance. Indeed, it can be easily verified 
that the above and below-threshold limits, Eqs. (jS]) and 
PTjl . of the pair resonance are reproduced by Eq. (fT^ . 
For V ^ Vc, Eq. ([T^ constitutes an approximate inter- 
polation formula, due to the attachment of an energy- 
independent width r to the two-particle resonance. 
We compare Eq. (|13p with a single-particle resonance 
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2e 



de 



/(e~f )-/(6+f ) 
{e-eaY + {T/2f 



(14) 

The qualitative difference between the two expressions 
arises from the appearance of the Fermi-liquid phase 
space factor [(ttT)^ + ~ cT^)^] in the two-particle res- 
onance Eq. This phase space factor can be re- 
moved by taking two derivatives with respect to voltage 
of Eq. p4|) . In this way, we find the relation 



rL) d xiPB.^2V,e,~eV/2], (15) 



dv^ 



(27r)2e2T/4 dV 



with the explicit replacements V — > 2V and ed — > Cd — 
eV/2. In view of the known properties of the single- 
particle resonance, this result constitutes our principal 
prediction. For T ^ T, Eq. p5|) predicts a Lorenzian 
peak in d^X/dV^ inside the Coulomb blockade diamond. 
Importantly, at ^ = K both d^J/dV^ and d^X^^^/dV^ 
have the same order of magnitude ~ T'^/{hVc). These 
results are illustrated in Fig. [H Note that, for T » 
r, Eq. (lisp also predicts temperature broadening of the 
peak in d^X/dV^. 
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FIG. 4: (a) The two-particle resonance induces a peak 
(dashed line) in (fiX/dV^ of width max(r, F), centered at 
V — Vc. The full- line curve shows that the two-particle reso- 
nance can be observed on top of the smoothly varying single- 
particle background, (b) The single-particle contribution to 
^3jiPfl^^y3 becomes singular at V ^ 2ed- 

Metallic dots. — Our predictions for transport via a sin- 
gle LS also extend to metallic islands with essentially 
zero level spacing. Transport through these islands can 
be modeled by the electric circuit shown in the inset of 
Fig. [21 The corresponding stability diagram includes a se- 
quence of Coulomb diamonds |9j, cf. Fig.[2l It is straight- 
forward to see that for metallic dots, the boundaries of 
two-electron ionization translate into a sequence of inner 
diamonds, as shown in Fig. [21 

Discussion and conclusion. — Previously it was be- 
lieved that in the course of cotunneling through a LS, 
electrons from the source arrive at the drain one by 
one. Here we demonstrated that there exists a well- 
pronounced, although more delicate, transport regime 
where two-electron processes contribute to the current. 
We emphasize that this regime is captured by the stan- 
dard Anderson Hamiltonian Eq. ([I]). 

Intriguingly, our reasoning is easily extended to 
regimes associated with A^-particle ionization of the LS 
{N > 2). These induce additional boundaries in the 
stability diagram Fig. [H at even lower voltages eV < 



eVc' 



{N) _ 



2ed/{2N — 1). A naive estimate of the corre- 



sponding near-threshold behavior of the current gives ~ 
0(^V - Fi'^>)[r2^-i(ey - eyi'^^)2(^-i)/(el4^^^)4(^-i)]. 
However, destructive interference between different se- 
quences of TV-electron transitions might lead to further 
reduction of the current. 

Throughout this paper, we considered an empty LS at 
zero bias (non-Kondo valley > 0). The analysis of 
the ionization process of the occupied LS (Kondo valley 
ed < 0) is entirely analogous, and differs only by the order 
of virtual transitions. 

This work was supported in part by the DFG through 
Sfb 658 and Spp 1243 (FvO), DIP (FvO and YO), 
ISF and BSF (YO) as well as KOSEF, KRF-2005-070- 
C00055, KRF-2006-331-C00118 (HSS). One of us (FvO) 
gratefully acknowledges hospitality by the Weizmann In- 
stitute, made possible by the EU - Transnational Access 
program (RITA-CT-2003-506095). 



5 



[1] M. A. Kastner, Rev. Mod. Phys. 64, 849 (1992). 
[2] L. Kouwenhoven and C. Marcus, Physics World 11, 35 
(1998). 

[3] L. Kouwenhoven and L.I. Glazman, Physics World 14, 33 
(2001). 

[4] J. Park et al, Nature (London) 417, 722 (2002); W. Liang 
et al, ibid. 417, 725 (2002); A.N. Pasupathy et al, Science 
306, 86 (2004). 



[5] M. Bockrath et al, Science 275, 1922 (1997). 
[6] A. Kamiuski and L.I. Glazmau, Phys. Rev. Lett. 86, 2400 
(2001). 

[7] P.W. Anderson, Phys. Rev. 124, 41 (1961). 

[8] A. Maassen van den Brink, G. Schon, and L. J. Geerlings, 
Phys. Rev. Lett. 67, 3030 (1991); P. Hadley et al, Phys. 
Rev. B 58, 15317 (1998). 

[9] G.-L. Ingoki and Yu. V. Nazarov, in Single Charge Tun- 
neling, edited by H. Grabert and M. Devoret, NATO ASI, 
Ser. B, Vol. 294 (Plenum Press, New York, 1992). 



